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One of the major goals of neutrino astronomy is to explore the otherwise unknown fluxes and 
interactions of ultrahigh energy neutrinos. The existing neutrino telescopes look at three types of 
events: particle showers, muons, and taus. In this paper we discuss the dependence of the event 
rates on the neutrino nucleon cross-sections as we scale the cross sections, with energy, in different 
scenarios beyond the standard model. Our focus will be on the IceCube detector. 
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Introduction: Neutrino- nucleon cross sections at ultrahigh energies are unknown at present. These cross sections 
might deviate from their standard model (SM) values due to the turn on of any new-physics processes above a TeV 
center of mass energies. A number of neutrino telescopes are already taking data and others are under developmental 
stage. One of the major goals of these telescopes is to determine neutrino-nucleon cross sections above a PeV. As the 
f^i ■ neutrino fluxes are not known either, the cross section information is lost if one merely looks at the absolute event 
rates integrated over neutrino energy. To determine the cross sections one needs to calculate the ratios of the upward 
+-> ■ and downward event rates. The knowledge of the scaling of event rates with cross sections also helps differentiate 
among the cross section models. In this work we study the scaling of neutrino event rates with the cross sections. For 
\^ , details see Ref. Q . A number of authors have done similar studies 0, 0, S IS 01 • 

I ■ A simple analytic model: Here we give a very simple model for the dependence of neutrino event rates on the 
ff") \ neutrino-nucleon cross sections. We will discuss the case for two generic detectors: an underground volume detector, 
like IceCube and a planar detector placed horizontally on the surface of the Earth. The cross section dependence 
of the downward event rates is trivial. The downward shower rates scale as -jg- ~ a t , where at is the total cross 

■ section due to all the processes that produce showers. Similarly, the lepton event rates are given by j^- ~ cr cc , where 
7-H • cr cc is the charged current (CC) interaction cross-section. 

The dependence of the upward rates on the cross sections is derived as follows: i) the rates depend on the upward flux 
at the detector: ~ Fd, where Fd depends on the attenuation cross section. If one ignores both the down-scattered 
neutrinos and the coupling among different neutrino flavors due to tau decays, the flux at the detector at nadir angle 

9 is Fd = F e~ l / Xa . Here F is the flux falling on surface of the Earth; I is the length neutrino travels through the 
Earth to reach the detector, and A a = l/NApo- a (Na is Avogadro's number and p is the density of the Earth which 

' is assumed to be constant in our simple model). a a is the attenuation cross section defined by a a = X^J/jUi- The 
^H, sum is over different types of interactions that a neutrino undergoes during its propagation through the Earth; here 
O ' y is the average inelasticity of an interaction defined as: y — 1 for a v e or charged current interaction; for v T 
charged current interactions, although a neutrino is completely absorbed, it is regenerated by tau decays; hence the 
1 effective y for j/ t 's is smaller than 1. For neutral current interactions y is just the conventional average inelasticity of 
the interaction (y ~ 0.2 — 0.25 8]). Given these considerations, for standard model, a a (SM) = 0.7a t (SM) around 

10 PeV pj. ii) The rates depend on the interaction cross section: j^- ~ er^; for showers, Oi is the total cross section 
due to the processes that produce showers and for leptons it is the charged current interaction cross section, iii) The 

■ rates depend on the area projected perpendicular to the direction of incoming neutrino: j^- ~ A p . For a volume 
detector, A p is almost a constant; for a surface detector, A p is a function of the nadir angle; for example, for a planar 
detector sitting horizontally on the surface of Earth, the projected area is A p = Acos8, where A is the surface area 
of the detector. This difference in A p for a volume detector and a planar detector causes the rates, integrated over 
angle, to be different for the detectors as we discuss below. 

rr/2 

The above discussion implies -j^- ~ j- J dflA p Fd ~ / dQ sind A p e~ l I Xa which for a volume detector, placed on 

o 

the surface (depth d = 0) in our simple model, reduces to -j^- ~ 4^(1 — e~ 2R / Xa ), and for a planar detector placed 
horizontally on the surface of Earth it reduces to -j^- ~ ^-(1 — e ~ 2 - R / A ° — lQ e ~ 2 R/ x a.y Now wc can write the scaling 

dE^ ' ' A7 



expressions for the event rates in two different limits: For Ai <C 2R, we get: j^- ~ 4^ = cri/a a for a volume detector 
and j^- ~ 4^ = ^i/^t f° r the planar detector. For 3> 2_R, the Earth becomes transparent and both for volume 
and planar detectors we get ~ &i (same as downward event rates as expected). 

Detailed Simulation: We solve coupled Boltzmann transport equations for three neutrino flavors and the taus 0. 
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As an example, we use IceCube as our detector and use the WB model flux Fq — 6 x 10 -8 (GeV.cm~ 2 .s^ 1 ) / (E(GeV)) 2 
[Tof. We include tau decay and the down scattering terms in the Boltzmann transport equations. We count all showers 
produced within the detector volume and take the energy losses of muons and taus to be independent of their energies; 
we also ignore tagging of taus here; our arguments about the scaling of event rates are not sensitive to these details. 
We do the simulation assuming three cases for neutrino-nucleon cross section models as explained below. Figures 1 
and 2 show the results obtained from detailed simulation. For details see Ref. 111. 
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FIG. 1: Expected event rates in IceCube vs. a (see text). Pairs of curves, from top to bottom, are for 10 6 ' 5 < E v < 10 9 GeV 
in five bins of width Alog 10 E v — 0.5. Fig 1(a) shower rates while the Fig. 1(b) shows lepton rates with d=1.9 km. 

Fig. 1 shows the results for the case where we assume the following cross-section models: a cc ~ aa cc (SM) and 
One = a<?Nc{SM) (these models are motivated by the fact that the QCD saturation effects might alter the scaling 
of neutrino-nucleon cross sections with energy at ultrahigh energies 11]). For this case, a a — (1 + a)a a {SM) and Oi 
[showers, leptons] = [(1 + a)at(SM), (1 + a)a cc {SM )]. Therefore, our simple analytic model implies that the rates 
-rE- ~ fj/ca are independent of a for these cross section models. Let us compare this prediction with the detailed 
simulation shown in Fig. 1. Fig. 1(a) shows results for the shower rates for a detector with IceCube characteristics. To 
illustrate the importance of detector depth d, we show results for d = 1.9 km (IceCube depth) and d = in Fig. 1(a). 
Fig. 1(b) shows lepton rates for d — 1.9 km case. We see in Fig. 1(a): a) For higher cross sections and energies, even 
a depth of 1.9 km makes a significant difference in the event rates as compared to the d = case. This is explained 
as follows: For d = 1.9 km case, even for horizontal neutrinos the distance traveled through the Earth is I ~ 220 km; 
hence the rate -j^- ~ - ie - 22 °( fcm )/ A a i s an upper limit on how fast the event rates can grow with the cross section 

for d = 1.9 km case. At higher energies and cross sections, A a < 220(fcm) (e.g. a(new) = 10a(SM), E v = 10 9 GeV 
=> A a ~ 50 km), hence the above upper limit on the rate is a decreasing function of a for our cross section models. 
On the other hand, for d = case, j^- ~ <Ti/cr a from our analytic model. These trends are obvious in Fig. 1(a). b) 
Another extreme limit is the one at low energies and cross sections (upper left corner in Fig. 1(a)) where both cases, 
d = and d = 1.9 km, are indistinguishable. This is because, in this limit, the Earth is becoming transparent to the 
neutrinos, c) For the case between the 'a' and 'b' extreme cases, our analytic model gives j^- ~ <Ji/a a ; however, 
in Fig. 1, we see there is a slight increase in the rates with a which is due to the details of the simulation that are 
absent from our analytic model (see the discussion of Fig. 2 below). 

Fig. 2 compares the scaling of event rates with the cross sections, for the energy bin 10 6 ' 5 < E v < 10 7 GeV, for 
different types of cross section models. The purpose of this figure is to illustrate the importance of inelasticity of the 
cross sections in the scaling arguments. All the results are for IceCube with d = 1.9 km. There are three triplets of 
curves: upper, middle, and lower; these triplets correspond to shower, muon, and tau rates, respectively. Each triplet 
has three curves (dashed, solid, and dotted) corresponding to three classes of models (labeled in the legend in Fig. 2) 
for the cross sections. We use the following cross section models in Fig. 2: For the dotted curves, our cross section 
model is the same as in Fig. 1. For the solid curves (motivation is black hole formation in low scale gravity models 
[l2|). in addition to the SM neutral current and charged current cross sections, we introduce a new cross section with 
inelasticity = 1 and we normalize it to the SM total cross section at a = 1. This implies a a = a a (SM) + aat(S M) and 
a [showers, leptons] = [(l+a)a t (SM), a cc (SM)} => ai/a a [showers, leptons] = [(l + a)/(0.67+a), 0.7/(0.67+a)]. For 
the dashed curves (motivation is graviton exchange in low scale gravity models.), in addition to the SM neutral current 
and charged current cross sections, we introduce a cross section equivalent to the SM neutral current cross section and 
we normalize it to the SM total cross section at a = 1. This implies a a = a a {SM) + a0.25o- t (SM) and di [showers, 
leptons] = [(l + a)a t (SM), a cc (SM)} aja a [showers, leptons] = [(1 + a)/(0.67 + 0.25a), 0.7/(0.67 + 0.25a)]. 
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FIG. 2: Comparison of rates, for the bin 10 6 ' 5 < E v < 
10 7 GeV, for different cross section models. There 
are three triplets of curves: Top: showers. Middle: 
muons. Bottom: taus. 



Fig. 2 is in good agreement with our simple analytic model; the event rates scale like ai/a a . The slight deviations 
from the expected scaling law ( -j^- ~ <x;/er a ) are due to a couple of details that we have ignored in our simple 
model but have included in the simulation (unlike the simulation, our simple model assumes constant values for Earth 
density and cross section inelasticity, and it also ignores the down scattering of neutrinos). 

In the above study, the normalization of the event rates was not our concern. To get an estimate of the event rates 
in the SM for IceCube, we have done a more realistic simulation which includes energy dependent energy losses for 
muons and taus; for tau energy losses we use approximation III from Ref. |13| and we use energy dependent effective 
area and effective volume for IceCube as given in Ref. [_|. For WB model, we get the upward rates per year in the 
energy bin 10 65 < E„ < 10 7 GeV: showers = 0.7, muons — 0.8, and taus = 0.2 (without tag); tau rate with tagging 
(as defined in Ref. 9] ) is 0.1. The rates integrated over energy above a PeV are: showers — 2.4, muons — 2.2, 
taus = 0.9 (without tag), and taus with tag = 0.3. 

Summary: We have shown that, in general, the cross section dependence of the event rates depends on the event 
type, neutrino energy, and the cross-section model, a) Downward shower rate is -jS- ~ at] downward lepton rate is 

j§- ~ fee- b) At neutrino energies around 1-10 PeV and for a volume detector, within a decade of the change in 

cross sections around the SM value, upward shower rate is ^j— - ~ at/ a a and the upward lepton rate is ~ a cc / a a . 

For a surface detector, the upward lepton rate is ~ a cc / a\. c) For energies around a PeV and cross sections 
around 0.1 of the SM value, the Earth becomes transparent and the upward rates scale with the cross-section like 
the downward rates, d) As one goes to higher energies and/or cross sections, the detector depth becomes important; 
even a depth of d = 2 km can give a significantly different scaling of the rates with the cross sections, as compared to 
the d = case (Fig. 1(a)). e) Scaling of the rates with the cross sections is significantly different for different models 
of cross sections (Figs. 2). 
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